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We study directed rigidity percolation (equivalent to 
directed bootstrap percolation) on three different lattices: 
square, triangular, and augmented triangular. The first two 
of these display a first-order transition at p — 1, while the 
augmented triangular lattice shows a continuous transition at 
a non-trivial p c . On the augmented triangular lattice we find, 
by extensive numerical simulation, that the directed rigidity 
percolation transition belongs to the same universality class 
as directed percolation. The same conclusion is reached by 
studying its surface critical behavior, i.e. the spreading of 
rigidity from finite clusters close to a non-rigid wall. Near the 
discontinuous transition at p = 1 on the triangular lattice, 
we are able to calculate the finite-size behavior of the den- 
sity of rigid sites analytically. Our results are confirmed by 
numerical simulation. 



I. INTRODUCTION 

Central-force rigidity percolation (RP) p|-pl| is the 
mechanical equivalent of the usual percolation prob- 
lem ]22|-p6|]. In RP forces (vectors) must be trans- 
mitted instead of scalars. This problem has received 
increased attention recentl y, following the development 
of mean-field theories [ |l^|Jl^Jlq| as well as of power- 
ful combinatorial algorithms ]27| |30| for its numerical 
study [p~3| — p~7]j3~S| — - As a result of these efforts, a deeper 
understanding of the rigidity transition has emerged, al- 
though some open questions remain. 

Bethe lattice calculations for RP |ll|[l8| with an ad- 
justable number g of degrees of freedom at each site have 
been used to obtain the behavior of the spanning cluster 
density P oc (p) as a function of p, the dilution (bond or 
site) parameter. For g = 1 one has usual (scalar) perco- 
lation, displaying a continuous transition with f3 = 1. 
But for any g > 1, the order parameter has a dis- 
continuity at a finite critical value p c . Thus the rigidity 
transition is discontinuous for d — > oo Jl^Jl^]. Other 
MF approximations also predict a first-order RP transi- 
tion [ [I4 ] 

On triangular lattices on the other hand, there is a di- 
vergent correlation length and the RP transition is second 
order 13-[l7],|2(],|2l|], but in a different universality class 
than usual percolation [ pOpl| ]. Some of the numerical 
evidence in 2d is consistent with a small discontinuity in 
the order parameter P^ , or a very small value for /?, but 
the precise interpretation of this evidence is still a matter 
of debate |l7j . In three dimensions the rigidity transition 



is undoubtedly second-order |31j . It is at present unclear 
in which fashion the RP transition becomes discontinu- 
ous as the dimensionality increases. Is there something 
like an upper critical dimension for RP, beyond which it 
is first-order? Or does it get increasingly "first-order" 
(i.e. (3 — > 0) as d — * 00? This analysis is further com- 
plicated by the fact that the character of the transition 
is lattice dependent. Hypercubic lattices in which sites 
have d degrees of freedom each cannot be rigid if they are 
diluted, but are rigid if undiluted and if they have appro- 
priate boundary conditions. Thus on hypercubic lattices 
the RP transition is "trivially first-order" at p c = 1, in 
any dimension. 

Similar considerations apply to directed lattices. Bethe 
lattices are directed by construction, since there is only 
one path between any two given sites. On directed lat- 
tices, rigid connectivity takes a particularly simple form. 
Imagine a rigid boundary to which a site with g degrees of 
freedom must be rigidly attached by means of rotatable 
springs (central forces). Each spring, or bond, restricts 
one degree of freedom. Thus the minimum number of 
bonds required to completely fix this site is g. Propa- 
gation of rigidity on directed lattices is then defined in 
the following terms: a site (with g degrees of freedom) at 
"time" t is rigidly connected to a boundary at t = if and 
only if it has g or more neighbors at earlier times who 
in turn are rigidly connected to the boundary. Thus, in 
contrast to undirected rigidity, which requires complex al- 
gorithms |27] |3^] that presently limit the maximum sizes 
to approximately 1.6 x 10 7 sites |^l| , po| , directed rigidity 
percolation (DRP) can be studied by means of a simple 
numerical procedure, and on much larger systems. 

It is interesting to notice that, on any directed lat- 
tice, rigidity percolation is equivalent to bootstrap perco- 
lation (BP), a modified percolation problem in which a 
site bel ong s to a cluster if at least m of its neighbors also 
do [p2]-|40|. Bootstrap percolation on undirected lattices 
attempts to describe certain systems in which atoms be- 
have magnetically only if they are surrounded by a "large 
enough" number of magnetic neighbors. A second reason 
for interest in BP is the search for novel critical behaviors 
in percolation ]37| , |38[ |, but the present understanding of 
this problem indicates that BP is either "trivially first- 
order" with p c ~ 1 or second-order and in the universality 
class of scalar percolation |fS9|-[tl|. 

Early studies of semi-directed m = 2 BP on square 
lattices seemed to indicate a transition at a non-trivial 
p 1 34 3(J, but rigorous arguments J35|] later showed that 
Pc = 1 in this case. To our knowledge there are no pub- 
lished studies of directed bootstrap percolation (DBP) 
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displaying a second-order transition. 

It has been recently conjectured E3| that any continu- 
ous transition in a nonequilibrium process with a scalar 
order-parameter, and a non-fluctuating, non-degenerate 
absorbing state must be in the same universality class 
as directed percolation(DP) p2|-|5l|. According to this, 
DRP-DBP would belong to the DP universality class in 
all dimensions for which it has a continuous transition. 

It is thus interesting to study RP on finite-dimensional 
directed lattices of increasing dimensionality, both to test 
this conjecture and to understand in which fashion DRP 
becomes discontinuous as d — > oo. In this article we re- 
port on our results for directed rigidity percolation (DRP, 
equivalent to DBP) on several (l+l)-dimcnsional lattices 
displaying first- and second-order phase transitions. 
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w(t) ~ t* (2) 

where \ = v ^l v \\ an( l ^J- i s the critical index associated 
with the decay length £j_ of perpendicular or "space" 
correlations: £j_ ~ \p — Pc\~ v± ■ Averages are taken only 
over clusters still alive at time t. Finally, the average 
mass of a cluster grown from a finite seed at p c behaves 
as 

M seed {t)~t\ (3) 

where fj = {u\\ + Vy - f3 dens )/v\\ p5|. 

For comparison we also simulate numerically usual di- 
rected percolation (DP, which corresponds to g — 1). In 
the DP case a simple argument shows that f3 dens = p seed 
because of time-reversal symmetry: consider for simplic- 
ity bond dilution and choose an arbitrary point x at 
time t. Any configuration of occupied bonds connecting 
(x, t) to the boundary at t — 0, and thus contributing to 
P(p,t), when reflected in the time direction, contributes 
to P^ eed (t,p) if now a point-like seed is located at x. 
Since both the original and the time-inverted configura- 
tion have the same probability, P(p,t) = P^ eed (p,t) ex- 
actly for bond-diluted DP and therefore (3 dens = (3 seed . 
Notice that this equality implies fj + S — x = !• 



FIG. 1. Lattices used for directed rigidity percolation 
studies in this work, a) square lattice, b) triangular lattice, 
and c) five-neighbor lattice. All examples are shown here 
undiluted (p — 1). Rigidity propagates upwards in "time" 
from the rigid boundary (gray) at t = 0. A site is defined to 
be rigidly connected to this boundary if it has two or more 
neighbors at earlier times which in turn are connected. 

If the DRP transition is second-order, we pay particu- 
lar attention to the determination of critical indices as- 
sociated with the spreading of rigidity, as described in 
the following fid| , f45|| . As usual in the study of directed 
processes, we define D(p) to be the asymptotic density 
of "active" (rigid) sites, which is equivalent to the prob- 
ability P{p) that, at large times t, a randomly chosen 
point be rigidly connected to a totally rigid boundary at 
t = 0. If the dilution parameter p is lower than a critical 
value pc, rigidity does not propagate and P(p) = 0. If 
the transition is second-order, immediately above p c one 
has P(p) ~ (p - p c f " nB . 

If the evolution starts from a finite rigid cluster or 
"seed" at t = instead of a rigid boundary, one defines 
P£ eed (t,p) as the probability that the cluster grown from 
this seed still be "active" at time t. If the transition is 
second-order, P^ eed (t,p) ~ (p — Pc) 13 for p — > p+ and 
t — > co. At p c this quantity decays as 

pr d (t lPc )~t- s , (i) 

with S = f3 seed /z/y and v\\ the temporal (or parallel) cor- 
relation length exponent: £|| ~ |p~Pc| _1/|1 

The typical width w of a cluster grown from a finite 
seed at p c behaves as 
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FIG. 2. Propagation of rigidity from finite clusters. The 
starting configuration is a finite sequence of contiguous rigid 
sites (black dots), a) On the square lattice, the size of this 
rigid cluster shrinks in time, b) On the triangular lattice, the 
rigid cluster's size remains constant in time if the lattice is 
undiluted, but shrinks in time for any p < 1, thus p c = 1 for 
this case, c) On the undiluted 5n-lattice, rigid clusters grow 
in time. Therefore for this lattice there is a nontrivial value 
p c (depending on the type of dilution, i.e. bond or site) above 
which rigidity propagates forever. 

Although no such time-reversal symmetry exists for 
DRP, we find that f3 dens = f3 seed also in this case. Fur- 
thermore we find that DRP belongs to the DP universal- 
ity class, i.e. has exactly the same critical indices. Thus, 
there is no separate universality class for directed rigidity 
percolation as there is for undirected rigidit y p ercolation. 
This is consistent with a recent conjecture |42j according 
to which any nonequilibrium process with a single ab- 
sorbing state will belong to the same universality class 
as DP. 

We also stud- 

ied the surface critical behavior |g|-^j5§|| of DRP, 
by means of simulations in the presence of an absorbing 
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boundary. In the DP case, the presence of the absorbing 
wall is known to only modify the survival exponent (3 seed . 
We find that this is also the case for DRP, and the new 
exponent is also consistent with the one obtained for DP 
with a wall. 

In section |l| we present our numerical results for DRP 
on directed lattices, with and without absorbing walls, 
and estimate the relevant critical indices associated to 



fects are different since now propagation of rigidity can 
exist without periodic boundary conditions. We discuss 



this case in detail in Sec. Ill 



the second-order transitions. Section III describes DRP 
on a directed triangular lattice. This case has a first- 
order transition at p = 1, and can be solved exactly for 
p-»l. 



II. NUMERICAL SIMULATIONS 

In order to simulate DRP we store a binary variable 
per site, indicating whether the given site is or not rigidly 
connected to the boundary at t = 0. We use the by 
now standard techniques of multispin-coding (MSC) [ p2[ , 
which allow us to store 64 binary variables in an inte- 
ger word, and also to update all of them simultaneously. 
Since interactions are short-ranged in the time direction, 
we only need to keep in memory a maximum of three 
consecutive lines of the system and we do this by means 
of three linear arrays which are reused periodically. 

We have considered three different oriented lattices: 
square, triangular and 5n-lattice (see later). The first two 
of them display trivial behavior (i.e. rigid only at p = 1) 
and the third one presents a continuous DRP transition 
at a non-trivial p c . In the first place we discuss DRP on 
a square lattice as depicted in Fig. |p,. Each site at time 
t has two neighbors at time t — 1. This is the minimum 
number of neighbors needed for rigidity with g = 2 and 
thus any amount of dilution is enough to impede propa- 
gation of rigidity. Therefore square lattices are not rigid 
at any p < 1. If p = 1, rigidity propagates only if bound- 
ary conditions are appropriate (e.g. rigid, or periodic, 
but not open). Any finite rigid cluster of size / shrinks to 
zero in I time-steps, as shown in Fig. |^a. The same would 
happen on directed ci-dimensional hypercubic lattices if 
g = d. 

We next consider the triangular lattice, oriented as 
shown in Fig. [T)b. Each site has three neighbors at earlier 
times. Despite the number of neighbors being larger than 
the minimum required (two), this lattice is also unable 
to propagate rigidity if diluted by any amount. To see 
why this is so, consider Fig. 0b, where one starts from 
a finite cluster of rigid sites (black sites) at t = 0. If 
the lattice is undiluted (p = 1), this cluster would just 
propagate unchanged in "time" . If the lattice is diluted 
by any amount, this rigid cluster would gradually shrink 
and eventually disappear. Thus for this lattice p c = 1, 
the same as for the square lattice. In contrast, finite-size 
effects are expected to be quite strong on the triangular 
lattice, since the lifetime of a finite rigid cluster diverges 
as p — y 1, no matter its original size. Also boundary ef- 



In order to have a nontrivial p c for DRP, we use trian- 
gular lattices augmented with two further bonds per site. 
These extra bonds connect layers t and t — 2, as shown in 
Fig. |l|c. This makes a total of five neighbors per site and 
we call this the 5n-lattice for simplicity. Now consider 
what happens when starting from a finite cluster of rigid 
sites on an undiluted 5n-lattice. As shown in Fig. ||c, the 
size of the rigid cluster expands in time with a constant 
angle if p — 1. Thus there will be a nontrivial value p c , 
above which rigidity propagates forever. We find that 
the DRP transition is second-order on this lattice. For 
comparison we also simulate DP on the square lattice. 

We next discuss our numerical results for DRP on the 
5n-lattice, and compare them to DP on the square lattice. 
We typically start our simulations from a finite seed of 
contiguous rigid sites and let the system evolve for 10 5 
timesteps (or until all activity dies out) and measure the 
survival probability p^ eed : cluster width w and average 
mass M see( i as a function of time. 
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FIG. 3. Survival probability for DRP clusters grown from 
finite seeds on site-diluted 5n-lattices with (from bottom to 
top) p = 0.70400, 0.70480, 0.70490, 0.70500, 0.70505 (dashed), 
0.70510, 0.70520, 0.70550, 0.70600 and 0.70700. Averages were 
taken over 10 5 realizations on systems of width 3840 sites. 

In a first set of simulations we estimate the critical den- 
sity p c for DRP on site-diluted 5n-lattices, by measuring 
Pa £ed (t) at different values of p and identifying the one 
for which the asymptotic behavior is closest to a straight 
line in a log- log plot (Figure ||). From these data we es- 
timate p^ RP = 0.70505 ± 0.00005. In contrast to P° eed , 
which shows appreciable curvature for off-critical values 
of p, the slopes of the cluster mass M see d and the mean- 
dering width w{t) in a similar log- log graph show little 
variation when p =/= p c . For DP on site-diluted square 
lattices, we use the estimate M p c = 0.64470. 
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FIG. 4. Critical behavior obtained from finite-seeds at 
the second-order transition for DRP on site-diluted 5n-lattices 
and for DP on site-diluted square lattices, a Survival proba- 
bility, b Cluster mass, c Meandering width. 

Figure ||a shows P se ed (i) for DRP (5n-lattice) and DP 
(square lattice) at their respective critical values. As- 
suming power-law corrections to (Q), we fit P se ed(Pc, t) — 
at- s (l + bt~ u ) and find 6 DRP = 0.15 ± 0.01 and S DP = 
0.16 ±0.02. 

The cluster mass M(t) and the meandering width w(t) 
behave as shown in Figures 0b and c respectively. From 
these data we estimate fj DRP ^= 1.47±0.01, t) dp = 1.47 ± 
0.01, X DRP = 0-633 ± 0.005 and X DP = 0.631 ± 0.005. 

These estimates are consistent with the more precise 
values @,@ S DP = 0.1594, f] DP = 1.4732 and X DP = 
0.6327, suggesting that DRP and DP are in the same 
universality class. 
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FIG. 5. Same as Fig. but now finite clusters start near 
an absorbing wall at x = 0. a Survival probability, b Cluster 
mass, c Meandering width. 

In order to further test whether DRP has the same crit- 
ical behavior as DP, we also studied DRP in the presence 
of an absorbing wall (DRPW). For DP with an absorb- 
ing wall (DPW) it is known g|-^,|^,|l| that the sur- 
vival exponent p seed is replaced by (3f , while z/|| and 
v\_ remain unchanged. Therefore only 6 is expected to 
change due to the presence of the absorbing wall. Our 
results are displayed in Fig. [|, and from them we ob- 
tain 5 DRPW = 0.423 ± 0.003, fj DRPW = 1.48 ± 0.01 and 
x drpw = o.62 ±0.02. Notice that 5 DRPW , fi DRPW and 
— X DRPW no longer add up to one since (3 dens and /3f eed 
are independent exponents. These results are entirely 
consistent with the values obtained for DPW by other 
authors MM. 
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III. DRP ON THE TRIANGULAR LATTICE 

This case is marginal as already advanced, since any 
amount of dilution will destroy rigidity and thus p c = 1, 
but on the other hand the lifetime of finite clusters is not 
finite as on the square lattice, but diverges as p — ► 1. As 
we show now, it is possible to obtain finite-size effects 
analytically for q = (I — p) << 1. 




a b 



FIG. 6. The directed triangular lattice presents a 
first-order transition at p = 1, which can be described in sim- 
ple terms. Arrows indicate absent sites, assuming site dilu- 
tion. When p < 1 a) defects (pairs of adjacent non-rigid sites) 
nucleate at rate p(p) per unit length, giving rise to non-rigid 
regions that b) widen in time with velocity v(p). 

Assume one starts from a completely rigid boundary at 
t = 0, on a triangular lattice of infinite width (Fig. |6p,). 
Let q = (1 — p) « 1 be the dilution parameter. We do 
not need to specify for the moment whether we are deal- 
ing with bond or site dilution. For short times all sites are 
rigidly connected to the lower boundary, but soon some 
non-rigid sites, or "defects" , will appear in the presence of 
either bond or site dilution. The smallest possible defect 
is a single non-rigid site, which happens with probabil- 
ity q per site on site-diluted lattices (one missing site) 
and with probability 3q 2 p + q a rs 3q 2 per site on bond- 
diluted lattices (two or three missing bonds) . This single 
defect "heals" immediately since each site above this one 
has three predecessors, but needs only two rigid ones in 
order to be itself rigid. 

A non-healing defect (in the following simply a defect) 
is created if two sites connected by a diagonal bond are 
simultaneously non-rigid, as in Fig. ^]a. All sites directly 
above these will have only one rigid neighbor and thus 
fail to be rigid, creating a "defect wall". Assume that 
these paired defects are nucleated with density p(p) per 
unit length (we calculate p later), and consider now the 
time evolution of the resulting defect wall. 

In the absence of dilution (q — 0), the boundaries of a 
non-rigid region stay unchanged in time (Fig. |]a). Rigid 
sites directly on this boundary have only two bonds (the 
minimum required number since g — 2) to rigid sites at 
earlier times. If one of these boundary sites fails to be 
rigid, all sites above it will also not be rigid. In this case 
the rigid boundary is displaced by one unit, as shown 



in Fig. |6|b. Therefore, for small but nonzero q, the rigid 
wall in Fig. ^b moves rightwards with an average velocity 
v = dx/dt which equals the probability for a boundary 
site to fail to be rigid. 

Neglecting fluctuations, we have a picture in which de- 
fects appear at a rate p(p) per unit length, giving rise 
to non-rigid regions which widen in time with constant 
velocity v(p). The system will become completely non- 
rigid when all defect regions have coalesced, as depicted 
in Fig. 0. This picture of the rigid-non-rigid transition 
is related to the Polynuclear Growth model (PNG) |53j , 
which has been extensively studied in the area of crystal 
growth. For our discussion of DRP we only need a few 
results which can be derived by means of simple argu- 
ments. 
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FIG. 7. Defects (empty circles) nucleate at rate p(p) per 
unit length, and the resulting "defect walls" widen in time 
with velocity v(p), eventually covering the entire system in a 
time of order t*. 

Assuming one knows the cone angle v(p) and the defect 
density p(p), it is easy to calculate the density of rigid 
points P(p,t) after t timesteps on a system of infinite 
width. A point (x, t) will be rigidly connected to the 
rigid boundary located at t = if it has not suffered 
the effect of any defect. In other words, if no defect 
has nucleated inside a "cone" with downwards opening 
angle v and whose vertex sits at (x,t). Let O = vt 2 be 
the area of this cone. Since defects nucleate randomly 
in space-time with density p(p), their number inside any 
given area O is a Poisson-distributed random variable 
with average Op. Thus 

P(p,t) = e-^\ (4) 

where 

t* = M~ 1/2 (5) 

is a characteristic time for the disappearance of rigidity, 
on an infinitely wide system. Using similarly simple ar- 
guments it is easy to see that the mean lifetime of a finite 
rigid cluster diverges as w _1 as p — > 1. 
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IV. CONCLUSIONS 




q = (1-p) 

FIG. 8. Mean time for the disappearance of rigidity, on 
site-diluted triangular lattices of width w = 128 (circles), 256 
(squares) and 512 (diamonds), as a function of q = 1 — p. 

We now calculate v{p) and p{p), and compare the re- 
sulting prediction for t* with our numerical results. Un- 
der site-dilution, the probability per unit time for a rigid 
wall to be displaced by one unit is simply v S it e — q. If 
bonds are diluted instead, one gets Vbond — 1 — p 2 ~ 2q. 
In order to calculate p, we notice that a pair of contigu- 
ous absent sites appears with probability p s u e — 2q 2 per 
site and per unit time on site-diluted lattices. On bond- 
diluted lattices on the other hand, creating such a pair 
requires at least three missing bonds. Thus pbond ~ q 3 - 
Finally one has (Eq. (§)) t* ite ~ q~ 3 / 3 and t* bond ~ q~ 2 . 
Figure || shows t* as measured on site-diluted lattices. 
These values are obtained by integrating in time the den- 
sity of rigid sites P(p,t). For an intermediate range of 
q, it is found that t* ~ q~ 3 ' 2 as predicted for infinitely 
wide systems. 

On a system of a finite width w with periodic boundary 
conditions, a crossover to a width-dominated behavior is 
expected if w « t*v, equivalently if w << (v/p) 1 ^ 2 , 
whereupon 

P(p,t) finite =e-*/*We , (6) 

and 



We considered directed rigidity percolation (DRP) 
with two degrees of freedom per site on three different 
(l+l)-dimensional lattices. This problem is equivalent 
to directed bootstrap percolation (DBP) with m = 2. 
On the square lattice, the system is only rigid at long 
times if p = 1. On triangular lattices a similar situa- 
tion happens, but this case has a non-trivial behavior for 
p — > 1, which we calculate analytically and confirm by 
numerical simulation. The mean lifetime of rigidity on 
infinitely wide systems is found to diverge when p — > 1 
as (1 — p)~ 3 / 2 for site dilution, and as (1 — p)~ 2 for bond 
dilution. The mean lifetime of a finite rigid cluster di- 
verges on the other hand as (1 — p)^ 1 in both cases. 

By augmenting the triangular lattice with two fur- 
ther bonds we define the 5n-lattice, which has a con- 
tinuous transition at pf RP = 0.70505 ± 0.00005 for site- 
dilution. We measure the critical indices associated with 
the spreading of rigidity and find that the DRP transi- 
tion belongs to the directed percolation (DP) universal- 
ity class, as a recent conjecture would indicate. A similar 
numerical study of DRP with an absorbing wall gives ex- 
ponents equally consistent with those of DP. Thus, while 
(undirected) rigidity percolation does not belong to the 
same universality class as usual percolation, the introduc- 
tion of directedness makes these two problems essentially 
equivalent at their respective critical point, i.e. on large 
scales. On (d+1) directed lattices, Bethe lattice calcu- 
lations indicate that the DRP transition becomes first- 
order for large d, while DP is always second-order (with 
mean-field exponents above its upper critical dimension 
d c = 5). We are presently extending this study to larger 
values of d. 
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^finite = OP) 1 ( 7 ) 

This regime corresponds to the defect-free area fi be- 
coming essentially a rectangle of height t and width w 
instead of a triangle of height t and base vt. According 
to equation (^), one should expect t* ~ q~ 2 for small q. 

This regime is observed for w — 128, but is less clear 
for larger values of w. Observation of this crossover for 
wider systems is numerically difficult, since it requires 
one to simulate very small q values, which makes the 
mean rigid times too large. 
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